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Abstract 

Consider sequential packing of unit volume balls in a large cube, in any dimension and 
with Poisson input. We show after suitable rescaling that the spatial distribution of packed 
balls tends to that of a Gaussian field in the thermodynamic limit. The results cover related 
applied models, including ballistic deposition and spatial birth-growth models. 

1 Introduction 

The following prototype random packing model is known as the basic Random Sequential Adsorp- 
tion Model (RSA) for hard spheres on a continuum surface. Open balls -Bi, ra , i?2,n-.., of unit radius 
arrive sequentially and uniformly at random in the d-dimensional cube Q n having volume n and 
centered at the origin. Let the first ball B\^ n be packed, and recursively for i = 2, 3, . . . , let the i-th 
ball Bi yTl be packed iff Bi iU does not overlap any ball in Bi jU , ...,-Bi_i „ which has already been 
packed. If not packed, the i-th ball is discarded. Given a positive integer k, let N({Bi^ n , Bk t n}) 
be the number of balls packed out of the first k arrivals. N n ^(k) :— N{{B l n , B/~ n }) are called 
random packing numbers. 

Lattice packing is defined analogously to continuum packing, save for the obvious constraint 
that centers of incoming balls are constrained to lie on a lattice. 
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There is a vast literature involving versions of the RSA model on continuum and lattice sub- 
strates. There is a plethora of experimental results and a notable dearth of mathematically rigorous 
results, particularly in more than one dimension. For surveys of the extensive literature, see Evans 
[8], Senger et al. [23], Bartelt and Privman [4], Adamczyk et al. [1], Talbot et al. [26] and [21]. 

In addition to their fundamental role in adsorption modelling, sequential packing models arise 
in the study of hard core interactions in physical and materials science, spatial growth models 
in crystallography and biology (Evans [8], sect. Ill, Garcia-Ruiz et al. [9]), and in the study of 
polymer reactions [9, 21]. In modelling communication protocols (Coffman et al. [6]), RSA is 
called on-line packing. 

Consider input of size V(tti), that is the deposition intensity is a Poisson random variable with 
parameter r. Deposition intensity in the continuum is the average number of particles arriving 
per unit volume, whereas in the lattice it is the average number of arriving particles per lattice 
point. In [20], the authors show that for both continuum and lattice packing, the Poisson packing 
numbers N n ^{V(Tn)) satisfy a thermodynamic limit as well as a central limit theorem. 

Theorem 1.1 (LLN and CLT for packing numbers [20]) For all r e (0,oo) and all d > 1 there 
are positive constants ad, T and rjd t T such that 

N n , d (V(Tn)) 
n 

while 



ad, T c.m.c.c. (1.1) 



N n AV(rn))-EN nA V(rn)) ^ ^ ^ ^ 

and 

n-W^N^iVirn})^^. (1.3) 

Here c.m.c.c. denotes convergence of means and complete convergence, a<i,T is called the cov- 
erage function, and 7V(0, a 2 ) denotes a normal random variable with mean zero and variance a 2 . 
The limits (1.1-1.3) also hold if the Poisson packing numbers N n ^{V(Tn)) are replaced by the 
finite input packing numbers AT n ^([m]) [20]. 

Our goal here is to investigate the large n distribution of packed balls on the d-dimensional 
substrate Q n . Even in dimension d = 1 little is known about the distribution of the point process 
on R induced by the packed balls. Our main result is that the basic RSA packing process on 
Q n C R d induces a point measure which, when suitably rescaled, converges to that of a Gaussian 
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field on R d in the thermodynamic limit. The CLT given by (1.2) of Theorem 1.1 is a by-product. 
Our main result also holds for variants of the basic packing process. 

The packing process is defined on the infinite substrate R d in a natural way. In this context, 
we represent the centers of the incoming balls, together with their arrival times, as a point set in 
R d x R+. Points of R d x R+ are generically denoted by to := (x, t x ), where x eR d , t x e M+. Given 
a locally finite point set X C R d x R + , we let a(X) denote the subset of X which is accepted in 
the packing process. We let n(X) C R d denote the projection of a(X) onto R d . n(X) is the point 
process on the substrate R d formed by the accepted balls. 

Let V T denote a rate one homogeneous Poisson point process on R d x [0,r]. In the lattice 
setting, Q denotes a collection of rate one homogeneous Poisson point processes on R indexed by 
the lattice points Z d and embedded in a natural way into R d x R+ . V T and Q henceforth represent 
the input for the continuum and lattice packing processes, respectively. The processes n(P T ) and 
7r(Q) are thinned Poisson point processes; we will be interested in their spatial distribution. 

For any set A C R d , let V t .a be the restriction of V T to A x [0, r] and let Qa denote the 
restriction of Q to A. 

For any Borel set B C R d define the random fields 

v T (B):= ]T S X (B) 

and 

v T , A (B) := S *( B y> 

that is, v T is the random field on subsets of R d induced by the packing process on R d whereas v Ti a 
is the random field on subsets of A induced by the packing process on A. Since sup B v T {B) = oo 
and sup B v t ^a{B) < volume(A), we call v T and v t ^a the infinite volume and finite volume packing 
measures, respectively. Edge effects show that in general these two random fields do not coincide 
on subsets of A. 

We similarly define random fields induced by the lattice packing process as follows. For any 
Borel set Bel 11 define 

/i(B):= Yl 5 ^ B ) 

xe-rr(Q) 

and 
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Consider the rescaled infinite volume continuum packing measures 



; _ Ut(\B) -Eu r (\B) 



and the rescaled finite volume continuum packing measures 



Vt,a,\K b ) — • I 1 - 5 ) 



The rescaled centered lattice packing measures [i\ and fj,A,\ are defined analogously. All pack- 
ing measures fi (or v) considered here are spatially homogeneous in the sense that the vectors 
(u(Bi), v{Bk)) and {v(B\ + x), v[Bk + x)) have the same distribution for any x G K d and 
any Borel subsets B\, Bk- 

In what follows all random measures are defined on the Borel subsets of R d . Recall that 
a sequence of random fields /x n , n > 1, converges in distribution to /x if and only if all finite 
dimensional distributions ([i n (Bi), yu„ converge to (fi(Bi), ...,n(Bk)), where Bi, ...,Bk are 
Borel subsets of R d . The following are our main results. 

Theorem 1.2 (Infinite volume packing measures converge to a Gaussian field) 

(i) (Poisson input continuum packing) For all t < oo, v Tt \ converges in distribution as A — > oo 
to a generalized Gaussian random field with covariance kernel K V T concentrated on the diagonal, 
that is 

K v , T (x, y) := C v , T 5{x - y), 

where C V ^ T is a constant depending on r. 

(ii) (infinite input lattice packing) fi\ converges in distribution as A — > oo to a generalized 
Gaussian random field with covariance kernel concentrated on the diagonal, that is 

Kn{x,y) : = Cf,S(x - y), 

where is a constant. 

The next theorem clearly extends Theorem 1.1. To obtain Theorem 1.1 we simply let A be the 
unit cube centered at the origin of K d . 

Theorem 1.3 (Finite volume packing measures converge to a Gaussian field) Let A C M. d have 
piecewise smooth boundary. Then: 

(i) (Poisson input continuum packing) For allr < oo, v T ,A,\ converges in distribution as A — > oo 
to a generalized Gaussian random field with covariance kernel K V , T . 
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(ii) (infinite input lattice packing) ha,\ converges in distribution as A — > oo to a generalized 
Gaussian random field with covariance kernel . 

Remarks, (i) We will see in the sequel that 



Jv 

where r\ and ri are the one and two point correlation functions for the spatial point process of 
packed points on R d . There is a similar expression for C^. Neither K V T nor depend on the set 



(ii) By definition of weak convergence, Theorem 1.2 tells us that for any Borel sets B\, ...,B m 
in M. d the m- vector with entries given by the i/ t ,a measure of B\, B m tends to a Gaussian limit 
with covariance matrix C;/, T volume(_Bi n Bj). Theorem 1.3 makes an analogous statement for the 
m-vector with entries given by the v T ,A,x measure of B\, B m . 

(iii) Proving convergence to a Gaussian field for infinite input continuum packing remains an 
open problem. In dimension d = 1, Dvorctzky and Robbins [7] showed that the number of packed 
balls asymptotically converges to a normal random variable. Theorem 1.3(ii) adds to results of 
Penrose [18], who shows that the number of packed balls in the lattice setting satisfies a central 
limit theorem. 

Acknowledgements. The authors are grateful for conversations with Y. Suhov, whose questions 
helped inspire the investigations in this paper. 



There are a multitude of variants of the basic RSA packing model [8, 21]. The approach taken here 
shows that the spatial distribution of packed balls (particles) for many of these variants converges 
to a Gaussian field. We discuss these variants below. Formal details may be found in [20]. 

2.1 Random shapes and types 

The basic RSA model assumes that the incoming particle is a ball of unit radius. However, this 
assumption may be relaxed to allow random shapes and random sets in R d , a widely considered 
model [4]. More generally, arriving particles may have a random type or spin, not necessarily 
representing shape or size. In one dimension such a model is considered by Itoh and Shepp [12]. 




A. 



2 Related Models 
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Instead of considering the point process denned by the point set of accepted particles, we may 
consider the point process consisting of only those accepted particles of a particular kind or type. 
Theorems 1.2 and 1.3 hold for the measures induced by such point processes. 

2.2 Time dependent models 

The basic RSA packing model can be generalized to include the case in which a packed ball 
remains in place for a random period of time at the end of which it is removed, i.e., desorbs. This 
is a dynamic model, which among other things, describes the reversible deposition of particles on 
substrates [23, 26]. Assuming that the spatial locations and arrival times of particles are given by 
a space-time Poisson process of unit intensity on M. d x [0,t], the contribution of a particle to the 
point process of adsorbed points is determined not only by whether it is accepted or not, but also 
by whether, if accepted, it desorbs by time r. The point measures given by the process of adsorbed 
points are defined analogously to (1.4, 1.5) and in the large A limit satisfy the Gaussian structure 
results of Theorems 1.2 and 1.3. 

We can also extend the basic packing model to a generalized version of the classical birth- 
growth model on R d in which cells are formed at random locations Xi G R d at times t i} i = 1, 2, ... 
according to a unit intensity homogeneous space-time Poisson point process. When a new cell is 
formed, its center Xi is called its "seed" . Once a seed is born, the cell around it initially takes the 
form of a ball of possibly random radius. If the initial radii of cells are zero a.s., then the model 
is referred to as the Johnson-Mchl model (see Stoyan and Stoyan [25]). The ball immediately 
generates a cell by growing radially in all directions with a constant speed. Whenever one cell 
touches another, it stops growing in that direction. New seeds, and the cells around them, form 
only in the uncovered space in R d . The point measures given by the process of seeds are defined 
analogously to (1.4, 1.5) and in the large A limit satisfy the Gaussian structure results of Theorems 

1.2 and 1.3. This adds to results of Chiu and Quinc [5], who only study the restriction to a 
large window of the infinite stationary birth-growth model on R d . Furthermore, as in Penrose and 
Yukich [20] (Theorem 2.1(b)), they consider only the number of seeds generated and not their 
spatial distribution. 

2.3 Ballistic Deposition Models 

The standard random particle deposition model considers random size i.i.d. (d + l)-dimensional 
balls ( "particles" ) which rain down sequentially at random onto a d dimensional substrate of volume 
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n, and centered around the origin. 

When a particle arrives on the existing agglomeration of deposited balls, the particle may slip 
and roll over existing particles, undergoing displacements, stopping when it reaches a position of 
lower height (the surface relaxation (SR) model [2]). If a particle reaches the surface R d , it is 
irreversibly fixed on it; otherwise, the particle is removed from the system and the next sequenced 
particle is considered. The rolling process does not displace already deposited particles; there is 
no updating of existing particles. For d = 1, the model dates back to Solomon [24]. Senger et al. 
[23] describe the many experimental results. 

The accepted particles all lie on the substrate and are represented by points in R d . The position 
of an accepted particle is a translate of the original location in R d above which it originally comes 
in. The point process of accepted particles defines a point measure on R d . The rescaled measures 
arc defined analogously to (1.4) and (1.5). The methods described below may be easily modified 
to show that in the large A limit, these measures converge to a Gaussian random field. This thus 
generalizes the results of [20] (Theorem 2.2(b)), which only considers the distribution of the number 
of accepted particles in the substrate. 

3 Auxiliary Results 

Throughout we let W := V x K+, where V :=R d . Recall that points of W are generically denoted 
by w :— (x,t x ), where x G V, t x G R+. The spatial location of the incoming balls, together with 
their arrival times, are represented by points in W. The interaction range of arriving balls is 
just the common diameter of the balls, that is to say equals two. Throughout, || • || denotes the 
Euclidean norm on R d . 

3.1 Exponential Decay 

Let X C W be a locally finite point process. Following [20] we make X into the vertex set of an 
oriented graph by including an edge from w\ := (x\,t Xl ) to W2 ■— {x2,t X2 ) whenever ||a;i— x%\\ < 2 
and t Xl < t X2 . Given w £ X, let A out (w, X) be the set of points (forward cone) in X that can be 
reached from w by a directed path in this graph (along with w itself). Let A- m (w,X) be the set 
of points (backward cone) in X from which the point w can be reached by a directed path in this 
graph (along with w itself). Finally, consider the "causal cone" 

^out,in 

(w, X) := A out (w, X) U A in (w, X). 
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The next result justifies describing the sets A ont ^ n (w, V T ) and A out ^n(w, Q) as "cones". 

Lemma 3.1 Fix r and let X be either V T or Q. There exist positive constants 7 := 7(r), (3 :— 
(3(t), and p := p(r), such that the causal cone associated with w := (x,t x ) € X belongs to the set 

C R {w) :={(y,t y ) :\\x-y\\ < /3\t x — t y \ + R} 

with probability at least 1 — pexp(— 7-R). 

Proof. For the case that X is the finite input set V T , this result is just Lemma 4.2 of [20]. For 
the case that X is Q, we will use an "invasion percolation" argument similar to that in [20]. The 
following argument actually also holds if X is the infinite continuum input V. First, we discretizc 
the space V by tiling it by lattice cubes of sufficiently large size so that no ball can intersect 
non-adjacent cubes. More precisely, we choose a full rank lattice L C V (a sublattice of Z d in the 
lattice case) with the elementary cube C, V = U^Ci, where C\ := C + i and the size of the cube 
C is chosen as to satisfy C + Bc2C (a further restriction in the lattice case is that the boundary 
of C docs not contain any point with integer coordinates, to avoid ambiguity). 

For each cube C\, i £ L, consider the arrival times of the points whose centers belong to this 
cube; they form a family P(i), i € L, of independent Poisson point processes. 

Now, fix a point w = (x, t x ) and assume, without loss of generality, that x belongs to the cube 
centered at the origin. Let us restrict attention to the forward cone and, to simplify notation, 
assume that t x = 0, i.e, w = (x, 0). Call a path w := no, tti, ■ ■ ■ ,tt p (of length p := \tt\) a collection 
of points in L, such that for any two consecutive points TTk and Wk+i, the cubes C T(e and C 7Tk+1 
share at least one vertex. For a path tt of length p, call the increasing sequence of arrival times 
= t < • ■ • < t p admissible if t k e P(itk) f° r all < k < p. Note that any path n admits many 
admissible sequences. Given a path tt of length p, let T(ir) be the infimum of the terminal values 
tp. 

The relations above imply that if the forward cone of (x, 0) does not belong to Cr(w), then for 
some i there exists a path n from to i such that |i| > f3T(n) + R, that is, 

T(n)<(\i\-R)/P. 

The probability that such a path exists is majorized by the sum over all possible paths 7r, 
starting at the origin, of the probability 

P[T(n) < (|7r w | - R)/0\ < P[T(n) < {k\n\ R)/0\ 
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(we use here the obvious fact that the distance between the last point of a path of length p and 
the origin is at most kp, where k is a constant depending only on the lattice L and the dimension 
d). 

The last piece we need before starting calculations is the fact that for a path it of length p, 
T(n) is simply the sum of p i.i.d. exponential variables Y. Choose 9 large enough so that 
£cxp(-6>y) < 3-( d+2 ) and choose (3 large enough so that e ek /P < 3. Let a := (kp - R)/0. Then 

P(tt) := P[T(w) <a}< e 0a (E[ e - eY ]) p . 

By definition of a and the choice of 9 and (3 we have that the above is bounded by 

e efcp//3 3 - P (d+2) < 3 - p (d+i)_ 

Since the number of paths starting at of length p is bounded by 3 dp , the probability that the 
forward cone of (x, 0) does not belong to Cr(w) is bounded above by the sum of the terms P(ir) 
over all possible paths n starting at and of length at least R/k. This probability is thus bounded 

by 

J2 3 dp P[T(n) <a}< J2 3 dp 3- p{d+1 \ 

p>R/k p>R/k 

which gives the desired exponential decay in R. 

Since the case of the backward cone may be handled by using a "time-reversal" , this completes 
the proof of Lemma 3.1. !|* 

Let X C W be a point set. In keeping with the terminology of statistical mechanics, we write 
a(w, X) = 1, if the ball centered at w is accepted with respect to X and otherwise. A trivial but 
useful observation is that if A[ n (x,Xi) = A[ n (x,X 2 ), then a(x,X\) = a(x,X 2 ). 

3.2 Correlation Functions 

To establish convergence of random point measures to a Gaussian field, one may employ the method 
of moments [11], which depends heavily on the use of correlation functions. 

Let V denote V T or Q. Recall that a(V) is the subset of V which is accepted in the packing 
process. Let 1b(c((P)) = 1 if ct(P) n B is non-empty, and let 1b((x(V)) = otherwise. Given 
Wi, ...,Wk G W, the fc-point correlation functions [22] of the point process a(V) of accepted points 
are defined as 

, x r E v l (wi+ein) a(T)---l (Wk+£k a ) a(T) 
rfc(wi '-" < Wfc):= £ ^o llm e^o £d+ i... £d+V ' 

1 k 
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where Wi ^ Wj,l < i ^ j < k, f2 is the ball in W of unit radius and oj is its volume. To obtain the 
definition in the lattice case one replaces the unit balls by the intervals {O} x [—1/2, 1/2]. 

To clarify the nature and existence of these correlation functions we introduce for w\, w k € W 
the functions 



r k (wi, ... ,w k ) := E v 



The functions r k are the probabilities that all points {wi, ... , Wk} are packed with respect to an 
independent sample from V . We notice here that this probability is always positive (unless the 
balls centered at the points x\, . . . ,x k themselves intersect), regardless of the heights t Xi of the 
points. 

We want to show that r k arc continuous in the arguments Wi := (xi, t Xi ), 1 < i < k. We will 
show the continuity of r k in the setting of the continuum V = V T ; continuity in the lattice setting 
follows from the arguments below and is easier since it only involves showing continuity in the time 
coordinate. 



Since 



is bounded by 



\r k (wi, ... ,w k )- r k (w[, . . . ,w' k )\ 



fc-i 



^2 \r k (w[, . . . ,w'^ 1 ,w l , ...,w k ) - r k (w[, . . . ...,w k )\, (3.2) 

»=i 

it will suffice to show that r k is continuous in each of its k arguments. We will show 

|r fe (wi,... ,w k ) -r k (w[,... ,w k )\ < C\w\ - w[\. 

The other summands in (3.2) are bounded similarly. For all x G V, let B{x) denote the ball in 
V of unit radius centered at x. B(w), w € W, is defined similarly. Without loss of generality, 
assume that \xi — x[\ < 1, so that B(x\) and B(x 2 ) overlap. 

Consider the event n»=i < J ( w i7'P) (T i w i-,T : ') E[i=2 cr ( u, i!^')- If ^is event occurs, then either 
the ball B(xi) is packed and B(x[) is not packed or vice versa. This means that the oriented graphs 
on V U {w\, . . . , w k } and V U {w[, . . . , w k } are different. The only way that this can happen is 
if V satisfies at least one of the following two conditions (assume without loss of generality that 

%Xl — tx[ ) • 

(a) the cylinder set 

[B(x 1 )AB(x' 1 )] x 
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intersects at least one ball B(w),w E V, whereas the cylinder set 

[B(xi) n B(x' 1 )\ x [0,t Xl ] 

does not intersect any ball B(w), w e V, or 
(b) the strip 

[B(a;i)UB(si)] x [t Xl ,t x{ ] 
intersects at least one ball B(w), aieP. 
Condition (a) implies that 

(i) the cylinder centered at (x\ + x[)/2 with radius 2 — \xi — x[\/2 and height t Xl does not 
contain any point from V, whereas 

(ii) the cylinder centered at (xi + x[)/2 with radius 2 + \x\ — x[\/2 and height t x > contains a 
point from V . 

Conditions (i) and (ii) imply that there is a set having volume of order \x\ — x[\ and containing 
at least one point from V. Similarly, condition (b) implies that there is a set having volume of 
order \t Xl — t x > \ and containing at least one point from V. The probability that a subset of W of 
small volume contains at least one point from V is of the order of the volume. It follows that the 
events (a) and (b) have probability of order \x\ — x[ \ + \t xi — t x > \, showing that the functions Tk 
are continuous. 



Rewrite the correlation functions as 

- k k 

nWe,n(a(7>)) I IIW ei n(P) 



lim E-p 

ei,... ,e k ^0 



x 



nti We 4 n(P) 



Li=l i=l 

We discard the event that any of the small balls {wi + £if2}^ =1 contains more than one point of 
V as an event of lower order magnitude. The first factor in the product is just Tk averaged over 
the possible positions of these solitary points in the small balls u>i + eiSl. Continuity of Tk implies 
that, up to negligible terms, one can replace this factor by the value of Tk at (101, ... , Wk)- The 
second factor in the product tends, obviously, to 

k 

Uh( Wi ), 

i=i 

where h{w) = h(t) is just the density of the underlying Poisson point process V with respect to 
the volume form on W (in our situation, h(t) is simply the indicator function of the interval [0, r]). 
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Summarizing, given w\,...,w k € W, the correlation functions are given by the following explicit 
formula 

k 

r k {wi,... ,w k ) =r k (wi,... ,w k )Y[h(wi). (3.3) 

i=l 

If a(V T )(B) represents the number of points in the set B cM. d x M+ which get packed then its 
moments are expressed in terms of the one and two point correlation functions via 

E[a(V T ){B)] = [ n(w)dw (3.4) 

JB 

and 

E[a 2 (P T )(B)} = / r 2 (w 1 ,w 2 )dw 1 dw 2 + / r^w^dwx. 

JBxB JB 

Also, we have (section 14.4 of Stoyan and Stoyan [25]) 

Var[a(V T )(B)] = / r 2 (w 1: w 2 ) ~ ri(w 1 )r 2 (w 2 )dwidw 2 + / r 1 (wi)dwi. (3.5) 
JBxB Jb 

The following is a key definition. 

Definition 3.1 A family of functions {rj : V J — > K}!^! exponentially clusters if for some positive 
constants Ak^Ckj one has uniformly 

\r k +i{wi,... ,Wk,w[,... ,iOj)-r fe («;i,... , w k )ri(w[, . . . ,w{)\ < A ki iexp(-C k ,id(w,w')), 

where d(w, w') is the distance between the sets {u>i, . . . , w k } and {w[, . . . , w\}, that is the minimum 
of the pairwise distances \wi — w'j\. 

Exponential clustering is also known as weak exponential decrease of correlations [15], or simply 
exponential decay of correlations. 

Let v be a spatially homogeneous point process on M. d . Then v defines a spatially homogeneous 
point process vl on the lattice Z d via 

vl ■= ^2 v(Q x )5 x , 

where Q x is the unit cube centered at x. Thus v L are probability measures on R z<i . Exponential 
clustering of the correlation functions of v implies the exponential clustering of the correlation 
functions for v^. By following the methods of Malyshev [15], who restricts attention to probability 
measures on {— 1, l} z<i , one can show that the normalized lattice measures 

v L (\B)-Ev L {\B) 

VL,\{B) ■ 
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converge in distribution to a generalized Gaussian field [15] on the Borel sets of R d (see Iagolnitzer 
and Souillard [11], especially p. 576). The lattice measures i>l,\ approximate the normalized 
continuum measures 

v{\B) - Ev(XB) 
v x (B) .= 



in that the mean and variance of \v(XB) — v L {XB)\ are of order A d_1 (see Proposition 4.1 below). 
Thus the continuum measures v\ also converge in distribution to a generalized Gaussian field. (An 
alternative approach involves working directly with the cumulants and proving that exponential 
clustering implies exponential decay of cumulants, automatically yielding the CLT we are seeking; 
see [14], [3]). 

Combining with the expression for the variance (3.5), we obtain the following result, the con- 
tinuum analog of Malyshev's CLT [15] for Gibbsian random fields. 

Theorem 3.1 (Gaussian CLT) Let v be a spatially homogeneous point process such that its cor- 
relation functions exponentially cluster. Consider the rescaled centered measures 

u{XB) Ev{XB) 

vx[B) = — — • 

Then, as A — > oo, v\ converges in distribution to a generalized Gaussian random field with covari- 
ance kernel 

K(x,y) :=C6(x-y), 

C := / r 2 (0 , x) - ri(0)n(x)dx + ri(0). 
Jv 

In other words, for any Borel sets Si,... ,B mi as A — > oo, the vector (y\(B\), ...,v\(B m )) 
tends to a Gaussian limit with the covariance matrix C(vol(Bj n Bj)), 1 < i,j < m. 

3.3 Process of packed points 

Recall that V denotes either V T or Q. To show Theorems 1.2 and 1.3, we will show that the 
correlation functions r^, k > 1, of the point process a(V) cluster exponentially. This will imply 
that the correlation functions rjj of the point process cluster exponentially, which, by Theorem 
3.1, gives the desired result. 

To show exponential clustering of the r^, k > 1, we first establish that a('P) is localized near 
V. This localization is of course obvious in the off-lattice case, as the i-support of the process is the 
bounded interval [0, r]. We cannot remove this cut-off as the current way of proof needs a rather 



where 
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rapid decay of the correlation functions with t which is lacking in the off-lattice case. Indeed, we 
have some basic estimates on bounds for the correlation functions. For example, the decay of n is 
polynomial: in dimension d = 1 it follows directly from the Renyi formula that r\(x,t) ~ t~ 2 and 
in higher dimensions r\ can be also shown to be at least of order The next lemma shows 

that the correlation functions decay exponentially in t in the finite input off-lattice case as well as 
in the infinite input lattice setting. 

Proposition 3.1 The correlation functions r k , k > 1, of the point process a(P) decay exponen- 
tially with t's, i.e. for any k there exists positive constants A k and Ck such that for w\, ...,w k € W 
we have 

r fe (wi, . . . ,Wk) < ^4feCxp(-C fe max(tj)). 

i 

Proof. Let Wi — (xi,ti), 1 < i < k. Consider first the correlation functions for a(V T ). 
Proposition 3.1 is clearly satisfied as in this case the left hand side of the inequality vanishes for 
maxi<feii > t. 

Now consider the correlation functions for a(Q), that is the correlation functions for the point 
process of accepted points on the lattice. Notice that, obviously, 

7k{wi, ... ,Wk)< ri(wi) 

for any 1 < i < k. Further, f\(wi) is just the probability that the point w\ = (x, t) is packed with 
respect to Q. However, this implies that all points w s = (x,st),0 < s < 1, can be packed with 
respect to Q, and therefore, that none of the points w s , s < 1, are present in the sample Q. This 
has probability exp(— at), where a > is the intensity measure of the interval {(x, s), < s < 1}. 
Thus n = f i decays exponentially with t and Proposition 3.1 is proved in the lattice setting. 18 

Proposition 3.2 Let Ek,i ■— Ekj(wi, ...,Wk,w[, ...,w[;T') be the event that the backward cones of 
any pair of points Wi and w'j with respect to V do not intersect. Then 

\r k+ i(wi, ...,w k ,w[, - r k (wi, ...,w k )ri(w[, ...,w ; ')l < Ckd(P[E k l }) 1/2 . (3.6) 

Proof. Define for any event E measurable with respect to the sigma algebra generated by V 



r k (u>i, ... ,w k ) := E v 



It suffices to show that (3.6) holds with r k replaced by r k . Now \r k — r k \ equals 



E'p 



< A k (P[E c })^ 2 , (3.7) 
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by Cauchy-Schwarz and the boundedness of r k . 

Observe that by the definition of E kl i, the random variables 

k i 
Y[a( Wl ,V)l Ektl and J| V)l Ek 

i=i j=i 

are independent. 

Hence the difference rf_£'/ (toi, Wk, w[, W;) — M (toi, Wk)rf M (wi, wj) vanishes and 
using the estimate (3.7) we obtain (3.6). H 
The clustering of the correlation functions is captured in 

Proposition 3.3 The correlation functions r k , k > 1, of the point process a(V) cluster exponen- 
tially. 

Proof. We fix t and use Lemma 3.1 describing the localization of causal cones. Let V = V T \ 
the proof for V = Q is exactly the same. The constants j,p, and (3 are as in Lemma 3.1. Let 
Wi := (xi,U), 1 < i < k. Let w'j := (x'^tj), 1 < j < I. 

Consider two sets {w\, . . . , Wk} and {w[, . . . , w[} at distance d. That is d := min \wi — w'j\. 
We distinguish two cases. 

(a) All times t\, . . . ,tk,t'i, . . . ,t[ are less than In this case with probability at least 1 — 
(hp) exp(— 7<i/4), the causal cone for each point Wi, 1 < i < k, belongs to the set Cj/^Wi) which 
is a subset of the cylinder of radius d/2 centered at x%. The same is valid for w'p 1 < j < I. Let E 
denote the event for which this is true. As corresponding cylinders for points from different tuples 
do not intersect, we conclude that E implies Ek,i, where Ek y i is defined as in Proposition 3.2. Now 
apply Proposition 3.2. 

(b) If at least one of the times (say, ti) is larger than jk, then by Proposition 3.1, both 
Tk+i{w\, . . . u>k,w[, w ; ') and rk(wi, . . . , u>k) decay exponentially in t\ and therefore decay expo- 
nentially in d as well. This finishes the proof. I 

We want to show the Gaussian structure for the process tt(V) by applying Theorem 3.1. tt(V) is 
a translationally invariant discrete (meaning that the distance between any two points is uniformly 
bounded from below) point process. 

To this end, we have to check that the correlation functions cluster exponentially. This is easy, 
given Propositions 3.1 and 3.3. Indeed, given xi, Xk G V, the correlation functions for ir(V) are 
given as 

rl(xi, . . . ,x k ) := / r k {{x\,ti), . . . , (x k ,t k ))dti . . .dt k . 

Jti,... ,t fc >o 
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To prove the clustering inequality for a fc-tuple (xi, Xk) and an i-tuple (x[, x[) at a distance 
d we split the integration domain (R + ) fc x (R + ) z into two subdomains, namely [0, d] k+l and its 
complement. Since all U, 1 < i < k, and t'j, 1 < j < k, are less than d in the first domain, we 
apply Proposition 3.3, using the fact that the ^-distance between Wi and w'j is at least d and 
the polynomial bound d k+l on the volume of the integration domain. This gives a bound which is 
exponentially decaying with d. 

In the second subdomain we apply Proposition 3.1, resulting in the estimate from above of this 
second part as 




A v .k,i exp(-Cy, fc ,zs)dvol({max(ti, . . . ,t k ,t[,... ,t[) < s}), 



where Ay,kd '■= Ak + Ai + Ak+i and Cy,k,i '■= niin(Cfc, Ci.Ck+i)- This is obviously exponentially 
decaying with d and thus the correlation functions r£ cluster exponentially. 

Now applying Theorem 3.1 we obtain the Gaussian CLT for the rescaled measures v T ^\. This 
completes the proof of Theorem 1.2. : ft 



4 Proof of Theorem 1.3 

We now deduce Theorem 1.3 from Theorem 1.2. We have already noted that because of edge 
effects, the rescaled point measures v T .\ and v T .A,\ are in general not equal. However, one can 
estimate the difference between the point measures v Ty \ and v T ,A.\- We will do this by using the 
exponentially clustering of the two-point correlation function to upper bound the variance of the 
difference. 

More specifically, for an open set A C M. d , let tt(V)a '■= 7r (^') H A and define the point measure 
n \A : ~ ^xeir(P\A)—n(P)\A ^A lSi tne difference between the packing process on XA and the 
infinite packing process restricted to XA. Similarly define ir^ A := J2xeTr(V)\ A —K(V\ A ) 

We want to estimate the variance of the number of points defining the supports of the point 
measures tt1^ a and n^ A . If we could show that this variance is of order o(A d ), A — > oo, then 
it would follow that the random measures tt('Paa), after centering and rescaling by \/A^, have 
asymptotically the same distribution as the centered and rescaled random measures tt(T')\Ai that 
is, have a generalized Gaussian distribution. 

To prove these variance bounds, we again resort to representation of the moments as integrals 
of polynomials of the correlation functions. In our case we will need just the expressions (3.4, 3.5) 
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for the first two moments. 

For all k we let denote the correlation functions for wf A . Observe the following: 

(a) The correlation function r'^w) := P[w e tt(V\a)\w € V\a] decays exponentially with time 
t for the same reasons that n does. 

(b) The correlation functions cluster exponentially. The proof goes along the lines of the 
proof of Proposition 3.3, and in fact involves conditioning on the same event. 

In other words, Propositions 3.1 and 3.3, together with the constants there, are valid for the 
point processes tt^a- 

The crucial property is the next one, which ensures that the point processes irf A have support 
which localizes near the boundary of XA. 

(c) The causal cone A out ^ n (w, V\a) of any point w with respect to the process V\a is a subset 
of A out ^ m (w,T\A) with respect to V and coincides with it if A outi i n (w, V) does not intersect the 
boundary d(\A). 

Therefore, the correlation functions rf of the processes nf A decay exponentially with the 
distance of ^-projection to d(XA): 

rf(w) < A d cxp{-C d d(x, d(XAj)}, w := (x,t) 

for some positive constants Aq,Cq- 
Together these properties imply 

Proposition 4.1 Assume that the set A is piecewise smooth. Then the variance and the mean of 
the number of points in either of the processes irf A is of order 0(X d ^ 1 ). 



Proof Start with the mean, which by (3.4), is given by 




By properties (a) and (c) above, the correlation function rf decays exponentially (say, as D exp(— Cd)) 
with the distance d from d(XA) C V x {0} C W. Hence the integral in question is bounded above 

by 

D / aq>(-C's)dB x (s), 
Jo 

where B x (s) := vol({x e M. d : d(x,d(XA)) < s}) and where {x e M. d : d(x,d{XAj) < s} is the 
s-neighborhood around the set XA. Since the leading term for the volume of the s-neighborhood 
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around XA has the form sX (1 ~ 1 Ca, were Ca is a constant depending on the set A, the above integral 
has leading term in A which is polynomial of degree (d— 1). 

To estimate the variance, we just add the clustering condition and use the variance formula 
(3.5) to obtain 

var(7r^ j4 ) = / [7*2" W2) — r+ (wi)ri (w2)]dwidw2 + / r^(w)dw, 

JWxW JW 

where r+ and r\ are the correlation functions for the process n^ A (an analogous formula is also 
valid for var(7r^ 4 )). 

The second summand is just the mean, and the integrand in the first integral decays in W x W 
exponentially with the distance from d(XA) C W x W. Once again applying the integral estimates 
as above, we obtain var(7r^" A ) = 0(A d_1 ). I 



5 Gaussian Fields and Total Edge Length Functionals 

The above methods extend existing central limit theorems involving functionals of Euclidean 
point sets, including those in computational geometry, Euclidean combinatorial optimization, and 
Boolean models. Existing central limit theorems ([19] and references therein) show asymptotic 
normality of functionals such as total edge length, total number of components, and total number 
of vertices of a specified degree. 

These functionals are canonically associated with point measures on R d . Given a graph G on a 
locally finite point set X, we associate to the total edge length functional the point measure defined 
by giving each vertex x e X a weight equal to one half the length of the edges in G incident to x. 

For example, given X a point set, let ux ■= vnn(x) be the point measure associated with the 
total edge length of the nearest neighbors graph on X. Thus, if G(X) is the nearest neighbors 
graph on X and if E (a;; G{X)) denotes the edges incident to x, then vnn(x) is defined as 

v N N{x)(B) := g 51 N' 

xeBnX e£E (x;G(X)) 

where |e| denotes the length of the edge e G G. Point measures associated with the total edge 
length of the Voronoi tessellation, minimal spanning tree, and sphere of influence functionals, are 
defined analogously. 

Let vnn ■— fNN,v (respectively, vnn,a ■— ^nn,a,v) denote the random measures associated 
with the nearest neighbors graph on the Poisson point process V (respectively, TT\A). Analogously 
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to the rescaled packing measures (1.4) and (1.5), define the rescaled (infinite volume) measures: 

/ r>\ . vnn(XB) - Ev NN (\B) 

and the rescaled (finite volume) measures 

, m vnn.xa(^B) - Ev nn ,\a{XB) 
vnn,A,\(B) ■= 7= • (0.2) 



Consider now the correlation functions for the point measure vnn(x)- Define the correlation 
functions of the point measures vnn by: 

E v v NN (w 1 + eiQ) ■ ■ ■ vnn{wu +e k n)) 
n{Wl > ■ ■ ■ ' Wk) ■= e^o lim £ ^o efi...^ ' 

In contrast to the packing measures, the point measures vnn do not evolve with time r. Instead 
of considering their localization properties via causal cones in R d x [0, r], we only need to consider 
their localization properties in M. d . This is accomplished by considering the notion of stabilizing 
measures [20]. Given a point measure v, write ux when v has support X on R d . Say that v 
is stabilizing exponentially fast if there is a ball Br := Br(Q), centered at the origin, with an 
exponentially decaying radius R, that is, P[R > t] < exp(— Ct), such that the difference measures 
defined by 

vvr\B R uAuo{ ■ ) - vvr\B R uA{ ■ ) 

are invariant for all finite A C R d — -Br(O). In other words, changes in the environment outside 
Bji(0) do not change the values of the point measure v inside Br(Q). 

By straightforward modifications to the proof of Proposition 3.3, it follows that if point measures 
stabilize exponentially fast, then their correlation functions will cluster exponentially, and therefore, 
the rescaled point measures converge to a Gaussian field. 

Now we may show [20] that the point measure vnn stabilizes exponentially fast. We thus have 
the following convergence result. An analogous convergence result holds for the rescaled finite 
volume measures (5.2). 

Theorem 5.1 (Infinite volume nearest neighbor measures converge to a Gaussian field) The mea- 
sures vnn,\ converge in distribution as A — > oo to a generalized Gaussian random field with co- 
variance kernel K concentrated on the diagonal, that is 

K(x,y) :=C5{x-y), 

where C is a constant. 
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Moreover, if we consider Voronoi tessellations, sphere of influence graphs, or minimal spanning 
tree graphs on Poisson point sets, then the canonically associated point measures associated with 
the total edge length functional are exponentially stabilizing [20] and they thus converge to a 
generalized Gaussian random field. 
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